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Abstract. We report on XMM-Newton observations of the "Black Widow pulsar", PSR B 1957+20. The pulsar's X-ray emis- 
sion is non-thermal and best modeled with a single powerlaw spectrum of photon index 2.03^ ^. No coherent X-ray pulsations 
at the pulsar's spin-period could be detected, though a strong binary-phase dependence of the X-ray flux is observed for the 
first time. The data suggest that the majority of the pulsar's X-radiation is emitted from a small part of the binary orbit only. We 
identified this part as being near to where the radio eclipse takes place. This could mean that the X-rays from PSR B 1957+20 
are mostly due to intra-shock emission which is strongest when the pulsar wind interacts with the ablated material from the 
companion star. 
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1. Introduction 

Until now, more than 1700 rotation-powered radio pulsars are 
detected. Among them are about 10% which are millisecond 
pulsars (MSPs) (Manchester et al. 2005). They form a separate 
population. The majority of them resides in Globular Clusters 
| (c.f. Bogdanov et al. 2006). MSPs are presumed to have been 
spun up in a past accretion phase by mass and angular transfer 
from a binary companion (Alpar et al. 1982). Only about one 
third of them are seen to be solitary. It is believed that they lost 
their companion, e.g. in a violent supernova event. All MSPs 
possess very short spin periods of less than 20 ms and show a 
high spin stability with period derivatives in the range * 10~ 18 - 
10~ 21 . MSPs are generally very old neutron stars with spin- 
down ages t = P/2P of ~ 10 9 - 10 10 years and low surface 
magnetic fields in the range B oc ^(PP) ~ 10 8 - 10 10 G. 

At present, about 50% of all X-ray detected rotation- 
powered radio pulsars are MSPs (c.f. Bogdanov et al. 2006 
and references therein). Among them an extraordinarily rich 
astrophysics binary system which is formed by the millisecond 
pulsar PSR B 1957+20 and its 0.025 M Q low mass white dwarf 
companion (Fruchter, Stinebring, & Taylor 1988). The binary 
period of the system is 9. 16 hours. The spin period of the pulsar 
is 1.6 ms which is the third shortest among all known MSPs. 
Its period derivative of P = 1.69 x 10~ 20 s s _I implies a spin- 
down energy of E = 10 35 erg s" 1 , a characteristic spin-down 
age of > 2 x 10 9 years, and a dipole surface magnetic field 
of B ± = 1.4 x 10 s Gauss. Optical observations by Fruchter et 
al. (1988) and van Paradijs et al. (1988) revealed that the pulsar 
wind consisting of electromagnetic radiation and high-energy 
particles is ablating and evaporating its white dwarf compan- 
ion star. This rarely observed property gave the pulsar the name 



black widow pulsar. Interestingly, the radio emission from the 
pulsar is eclipsed for approximately 10% of each orbit by mate- 
rial expelled from the white dwarf companion. For a radio dis- 
persion measure inferred distance of 1 .5 kpc (Taylor & Cordes 
1993) the pulsar moves through the sky with a supersonic ve- 
locity of 220 km/sec. The interaction of a relativistic wind flow- 
ing away from the pulsar with the interstellar medium (ISM) 
produces an H a bow shock which was the first one seen around 
a "recycled" pulsar (Kulkarni & Hester 1988). 

In 1992 Kulkarni et al. published a contour map of the 
X-ray emission of PSR B 1957+20 which was derived from 
ROSAT PSPC observations. Although this ROSAT data were 
very sparse in statistics it led the authors to predict faint diffuse 
X-ray emission with constant surface brightness to be present 
along a cylindrical trail formed when the relativistic pulsar 
wind expands into pressure equilibrium with the interstellar 
medium behind the nebula. The much improved sensitivity of 
the Chandra and XMM-Newton observatories made it possi- 
ble to probe and investigate the structure and properties of this 
unique binary system in much higher detail than it was possible 
with ROSAT, ASCA or BeppoSAX. A narrow X-ray tail with 
the extent of 16 arcsec and the orientation to the north-east was 
detected from it in deep Chandra observations by Stappers et 
al. (2003). 

Searching in ROSAT data for a modulation of the pulsar's 
X-ray emission as a function of its orbital phase revealed a 
suggestive but insignificant increases in flux before (at phase 
<p ~ 0.17) and after (at phase <p ~ 0.4 - 0.5) the pulsar radio 
eclipse (0 = 0.25) (Kulkarni et al. 1992). Taking Chandra data 
into account revealed a hint that the lowest and highest fluxes 
are located during and immediately after the radio eclipse, re- 
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spectively. The statistical significance of this modulation ob- 
served by Chandra, though, is only 98% and thus prevents any 
final conclusion on it (Stappers et al. 2003). 

In this paper we report on XMM-Newton observations of 
PSR B 1957+20 and its white dwarf companion. The paper is 
organized as follows. In §2 we describe the observations and 
data analysis. §3 gives a summary and discussion. 



2. Observations and Data Analysis 

PSR B1957+20 was observed with XMM-Newton on October 
31, 2004 for a 30 ksec effective exposure. In this observation, 
the EPIC-MOS1 and MOS2 instruments were operated in full- 
frame mode using the thin filter to block optical stray light. 
The EPIC-PN detector was setup to operate in the fast timing 
mode. Because of the reduced spatial information provided by 
the PN in timing mode we use the MOS1/2 data for imaging 
and spectral analysis of the pulsar and its diffuse X-ray neb- 
ula while the PN data having a temporal resolution of 0.02956 
ms allowed us to search for X-ray pulsations from the pulsar. 
All the data were processed with the XMM-Newton Science 
Analysis Software (SAS) package (Version 6.5.0). Spatial and 
spectral analyses were restricted to the 0.3 - 10.0 keV energy 
band while for timing analysis events were selected for the en- 
ergy range 0.3 - 3.0 keV. 

2.1. Spatial Analysis 

Figure 1 shows the combined EPIC-MOS1/2 image of the PSR 
B 1957+20 system. The image was created with a binning fac- 
tor of 6 arcsec and by using an adaptive smoothing algorithm 
with a Gaussian kernel of cr < 4 pixels in order to better 
make visible faint diffuse emission. The extent of the diffuse 
emission with its orientation to the north-east is about 16 arc- 
sec which is consistent with the previous result derived from 
the Chandra Observations. However, the detailed structure of 
the X-ray emission from XMM-Newton can not be as clearly 
seen as from Chandra due to the 10 times wider Point Spread 
Function (PSF) of XMM-Newton. 

Inspecting the XMM-Newton MOS 1/2 image two faint fea- 
tures (denoted as A and B) which contribute to only about 
3% of the total X-ray flux are apparent. In order to investi- 
gate whether this faint features are associated with nearby stars 
we inspected the Digitized Sky Survey data (DSS) and the 
USNO-B1.0 Catalogue for possible sources. These catalogues 
which are limited down to 22 mag (Krongold, Dultzin-Hacyan 
& Marziani 2001) and 21 mag (Monet et al. 2003) do not re- 
veal possible counterparts. These features are not seen in the 
Chandra image though (Stappers et al. 2003). 

2.2. Spectral Analysis 

Combined EPIC-MOS1/2 data of PSR B1957+20 were ex- 
tracted from a circle of 30 arcsec radius centered at the pulsar 
position (RA (J2000) = 19 /! 59 m 36' 5 .77, Dec = 20°48'15".12). 
The selection region contains about 85% of all source counts. 
Background photons were selected from a source-free region 
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Fig.l. XMM-Newton MOS1/2 image of the PSR B1957+20 
-WD system with contour lines overlaid. The contour 
lines are at the levels of (5.1,6.1,9.4,16.3,32.4,69.0) x 
10~ 6 cts s~'arcsec~ 2 . The position of the pulsar is indicated. 



near to the pulsar position. Response files were derived by us- 
ing the XMM-Newton SAS tasks RMFGEN and ARFGEN. 

After subtracting background photons, in total 338 sources 
counts were available for a spectral analysis. The extracted 
spectra were binned with at least 30 source counts per bin. 
Assuming that the emission originates from the interaction of 
the pulsar wind with the ISM or with the stellar wind we ex- 
pect synchrotron radiation to be the emission mechanism of the 
detected X-rays. To test this hypothesis we fitted the spectrum 
with a power law model. Indeed, this model describes the ob- 
served spectrum with a reduced ;f 2 of 1.09 (for 8 D.O.F). The 
photon-index is found to be a = 2.03*n|i. The column absorp- 
tion Nh is 8.0 x 10 20 cmT 2 . For the normalization at 1 keV we 
find 1.5 ^'xlO -5 photons keV cm^sec -1 (1-cr confidence for 
1 parameter of interest). The spectrum (data and model) and the 
fit residuals are shown in Figure 2. 

The unabsorbed X-ray fluxed derived from the best fit- 
ting model parameters is f x = 8.35 x 10~ 14 erg s^ 1 cm' 2 
and f x = 7.87 x 10~ 14 erg s~ l cm' 2 in the 0.3 - 10 keV and 
0.1 - 2.4 keV energy band, respectively. The X-ray luminosi- 
ties in these energy bands - calculated for a pulsar distance of 
1.5 kpc - are L x (0.3 - lO.OkeV) = 2.24 x 10 31 erg s' 1 and 
L*(0.1 - 2.4keV) = 2.12 x 10 31 erg s~\ respectively. The con- 
version efficiency L x /E in the 0.1 - 2.4 keV band is found to 
be ~ 2.12 x 10~ 4 . 

In order to check whether the spectral emission character- 
istics changes for photons detected in a smaller compact region 
of 10 arcsec radius at the pulsar position we applied a spectral 
fit to this events only. The encircled energy within 10 arcsec 
is 60 %. In total, 186 sources counts were available for spec- 
tral fits. We did not find any significant change in the spectral 
parameters than reported above. 
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Fig. 2. Energy spectrum of PSR B 1957+20 obtained from the 
XMM-Newton MOS1/2 data. The plot shows the X-ray spec- 
trum fitted with an absorbed power law model (upper panel) 
and contribution to the^ 2 fit statistic (lower panel). 



2.3. Timing Analysis 

The EPIC-PN camera observed the pulsar in the fast timing 
mode. In this mode the spatial and spectral information from a 
64 x 199 CCD pixel array is condensed into a one dimensional 
64 x 1 pixel array (ID-image), i.e. the spatial information in 
Y-direction is lost due to the continuous read-out of the CCD. 
The complete photon flux (source plus DC emission from fore- 
ground or background sources located along the read-out di- 
rection) is accumulated and collapsed in the final ID-image, 
severely reducing the signal-to-noise ratio of pulsed emission 
and preventing the detection of weak X-ray pulsations from the 
target of interest. 

In order to search for X-ray pulsations from PSR 
B 1957 +20 we extracted 1639 counts from the CCD columns 
33 - 41 in which the pulsar got located. In order to increase 
the signal-to-noise ratio we restricted the analysis to the energy 
range 0.3 - 3.0 keV. Below and beyond this energy band the 
accumulative sky and instrument background noise exceeds the 
contribution from the pulsar itself (c.f. Becker & Aschenbach 
2002). However, still 80 % of the counts are estimated to be 
derived from the background. 

The photon arrival times were corrected to the solar sys- 
tem barycentre with the BARYCEN tool (version: 1.17.3, JPL 
DE200 Earth ephemeris) of the SAS package. As the pulsar is 
in a binary we corrected for the orbital motion of the pulsar by 
using the method of Blandford & Teukolsky (1976). 

As millisecond pulsars are known to be extremely stable 
clocks we used the pulsar ephemeris from the ATNF Catalogue, 
/ = 622.122030511927 Hz and / = -6.5221 x lO^sec" 2 (at 
MJD = 48196.0) to perform a period folding. By using the Z 2 
statistics (Buccheri et al. 1983) with the harmonics number (n) 
from one to ten no significant signal was detected at the radio 
spin period extrapolated for the epoch of the XMM-Newton 
observation. Restricting the period search to the various smaller 
energy bands did not change the result. A pulsed fraction upper 



limit of 9% (l-<x) is deduced by assuming a sinusoidal pulsed 
profile. 

Arons & Tavani (1993) predicted that depending on the 
flow speed and the degree of absorption and/or scattering by 
the companion wind the X-ray emission from PSR B 1957+20 
increases by up to a factor of 2.2 at the orbital phase before and 
after the radio eclipse. In order to test this prediction we cre- 
ated a lightcurve by binning all events in bins of 1 .5 ksec width. 
With an effective exposure time of about 30 ksec and an orbit 
period of 9. 16-h the XMM-Newton MOS 1/2 and PN data cover 
roughly 83 % and 92 % of one binary orbit, respectively. Table 
1 lists the first and the last photon arrival times recorded by the 
MOS 1/2 and EPIC-PN detectors and the corresponding pulsar 
binary orbital phase. The lightcurves resulting from the EPIC- 
PN and MOS 1/2 data are shown in Figure 3. For the lightcurve 
deduced from the PN-data (upper solid curve) it is clearly seen 
that before the radio eclipse, i.e. between the orbital phase 0.1 
and 0.25, the X-ray emission increases. The emission in the 
highest bin is about a factor of 3.0 stronger than at other orbital 
phase angles before and after the radio eclipse. This is in agree- 
ment with the predictions made by Arons & Tavani (1993). 

Although it is the first time that a significant X-ray flux 
modulation from PSR B1957+20 is observed, the flux increase 
near the phase of the radio eclipse is virtually only seen in the 
PN data. Owing to the short observation time of 30 ksec, which 
is less than the time of one full binary orbit, the orbit angles 
0.30 - 0.38 and 0.18 - 0.36 are not covered by EPIC-PN and 
MOS 1/2 data at all. It thus turns out that it was a great fortune 
that the EPIC-PN, more or less by chance, covered the orbital 
phase range of the radio eclipse and provided us evidence for 
the flux enhancement while the MOS 1/2 detectors were already 
switched off. 

3. Summary and Discussion 

Interaction between relativistic pulsar winds which carry away 
the rotational energy of pulsars and the surrounding medium is 
expected to create detectable X-ray emission. Indeed, there are 
about 30 pulsar wind nebulae (PWNe) currently detected in the 
X-ray band (e.g. Kaspi, Roberts & Harding 2006, Gaensler & 
Slane 2006, Kargaltsev & Pavlov 2006). However, these PWNe 
are all powered by young and powerful pulsars with spin-down 
energies of more than ~ 3.6 x 10 36 erg s~ l . Until now, only 
two MSPs are known to have X-ray nebulae. They are PSR 
B1957+20 (Stappers et al. 2003) and PSR J2124-3358 (Hui 
& Becker 2006). Both of them have tail-like structures behind 
the moving pulsars. These trails could be associated with the 
shocked relativistic wind confined by the ram pressure of the 
ambient ISM. 

The XMM-Newton data of PSR B1957+20 have provided 
observational evidence for a strong dependence of the pulsar's 
X-ray emission on its binary orbital phase. It is the first time 
that a significant X-ray flux modulation from PSR B1957+20 
is observed. The emission near to the radio eclipse is sup- 
posed to be beamed in a forward cone because the shocked 
fluid is accelerated by the pressure gradient as it flows around 
the eclipse region. Relativistic beaming would tend to give the 
maximum flux just before and after the pulsar eclipse which 
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Table 1. List of the first and the last photon arrival times and its corresponding orbital phase of PSR B 1957+20 for the MOS1/2 
and PN detectors. Arrival times are corrected to the solar system barycentre. The orbital phase is measured from the time of 
ascending node. 



Date set 


1st photon 


last photon 


Duration (s) 




Time (MJD) Orbital Phase 


Time (MJD) Orbital Phase 




MOS1 


53309.9688 0.3605 


53310.2830 0.1830 


27143.4 


MOS2 


53309.9666 0.3548 


53310.2846 0.1871 


27467.1 


PN 


53309.9791 0.3873 


53310.3301 0.3063 


30327.4 



MJD 

53309.95 2h 4h 6h 8h 9.16h 





N 

V10S1/2 
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Fig. 3. X-ray emission from PSR B 1957+20 within 0.3 - 3.0 
keV as function of the pulsar's orbital phase ((f)). We mapped 
one complete orbital period of this system at the starting point 
of MJD = 53309.95, i.e. <f> = 0.32. <p = 1.0 corresponds to the 
ascending node of the pulsar orbit. The upper curve was ob- 
tained from the XMM-Newton EPIC-PN (background level at 
66 cts/bin). The lower lightcurve is obtained from the MOS1/2 
data. The gray strip between the orbital angle 0.21 - 0.29 in- 
dicates the eclipse of the pulsar. Phase bins with zero counts 
correspond to phase angles not covered in the observation. 

for PSR B1957+20 is at orbital phase = 0.25. Arons & Tavani 
(1993) predicted that the immediately downstream flow veloc- 
ity of the shocked pulsar wind from the shock area along the 
line between the pulsar and the ablating companion is about 
c/3 and is even higher behind the relativistic shock when it 
passes around the companion. The tenuous relativistic plasma 
accelerates as it flows around the companion, possibly passing 
through a sonic transition to leave the binary system with the 
velocity larger than c/ V3. Due to Doppler boosting, the prob- 
able post-shock velocities of the relativistic wind then suggest 
X-ray emission variation around the orbital phase by a factor 
between 1.3 (v = c/3) and 2.2 (v = c/ V3). This numbers were 
estimated without considering the effect of absorption and/or 
scattering within the binary. On the contrary, the X-ray emis- 
sion at the eclipse may be reduced because of the obscuration 
of the shock by the companion. 

However, given the limited photon statistics of the XMM- 
Newton data it is not possible to investigate any binary-phase 
resolved imaging and spectral variation as a function of orbit 
phase or to determine the exact geometry of the peak emis- 



sion. Such analysis would allow us to investigate with higher 
accuracy than currently possible whether the X-ray emission 
from PSR B 1957+20 is present during all orbit angles or vir- 
tually only near to the radio eclipse while diffuse X-ray emis- 
sion from the PWN is present at any orbital angle. Indeed, the 
later is suggested by the current data and a confirmation would 
have a severe impact on our understanding of the pulsar's X-ray 
emission properties. 

As the present XMM observation covers barely one binary 
orbit, we stress that it can not be fully excluded that the increase 
in photon flux near the orbital angle of the radio eclipse is due 
to a single "burst like" event or that the peak flux emission 
varies from orbit to orbit or on longer time scales. Clearly, a 
repeated coverage of the binary orbit in a longer XMM-Newton 
observation and a comparison with the 2004 data would answer 
this question immediately. This, in addition, would provide us 
not only a better photon statistic but would also allow us to 
determine the emission geometry with a much higher accuracy 
than currently possible. 
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